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Abstract 

We investigate Nuclear Lattice Effective Field Theory for the two-body system for several lattice spacings at lowest order in the 
pionless as well as in the pionful theory. We discuss issues of regularizations and predictions for the effective range expansion. 
In the pionless case, a simple Gaussian smearing allows to demonstrate lattice spacing independence over a wide range of lattice 
spacings. We show that regularization methods known from the continuum formulation are necessary as well as feasible for the 
pionful approach. 
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1. Introduction 

The Nuclear Lattice Effective Field Theory (EFT) method 
01 has led to impressive progress in the last decade and it has 
been applied to few- and many-body-systems successfully, for 
reviews see e.g. Refs. HI- The lattice spacing serves as a 
natural UV regulator for the theory, as for a given value of a 
the maximal momentum is /? max = n/a. Although these cal¬ 
culations give a quite good description for the phase shifts, en¬ 
ergy levels, etc., almost all calculations have been done for a 
fixed lattice spacing a ^ 2fm, corresponding to a soft mo¬ 
mentum cut-off of about 300 MeV. This allows one to treat all 
corrections beyond leading order (LO) in perturbation theory. 
However, the cut-off dependence or lattice spacing dependence 
has not been analyzed systematically and there are still some 
problems in the two-nucleon system like the relatively poor de¬ 
scription of the 3 Si- 3 Di mixing angle 01 . Further, such soft 
potentials seem to lead to some overbinding in medium-mass 
nuclei, as discussed in Ref. 0. Also, it has been shown that 
the leading order four-nucleon contact interactions need to be 
smeared to avoid a cluster instability when four nucleons reside 
on one lattice site 0. One might argue that the extension of 
such smearing methods also to the pion exchange contributions 
leads to a natural regularization of the lattice EFT, allowing to 
vary the lattice spacing freely but using an explicit momentum 
cut-off in the spirit of the work of Ref. |7[]. More precisely, this 
inherent physical cut-off was implemented by formulating the 
lattice action in terms of blocked fields. 

In this paper, we will focus on the neutron-proton two-body 
system at lowest order and discuss the lattice spacing depen¬ 
dence systematically. In addition, we discuss the necessity 
of regularizing the one-pion-exchange potential and provide a 
method that goes beyond smearing and is borrowed from con¬ 
tinuum calculations, which leads to the lattice spacing indepen¬ 


dence of observables for a broad range in a, see Ref. 0. 

While most of the calculations solve the transfer matrix using 
Monte Carlo methods or the Lanczos method for small eigen¬ 
values of large sparse matrices, we use here the Hamiltonian 
formalism and solve it with the Lanczos method. Using this 
approach we can eliminate the discretization in the time direc¬ 
tion and we have to consider only the variation in the position 
space discretization. In the following, all expressions are given 
in lattice units and one has to multiply the lattice results by the 
appropriate power of the lattice spacing a to get the physical 
values. Note also that we show simulations for various large 
enough volumes so that Liischer’s finite volume formulas are 
sufficient for the infinite volume extraction and we can entirely 
focus on the remaining dependence on the lattice spacing. 

In what follows, we will first display the necessary formal¬ 
ism to calculate the neutron-proton system to lowest order on 
the lattice. It is important to already improve the free Hamilto¬ 
nian so as to be as close as possible to the free non-relativistic 
dispersion relation. At very low energies, one can consider the 
theory with contact interactions only, the so-called pionless the¬ 
ory. As we will show, the smearing of the contact interactions 
can be used as a regulator, leading to regulator-independent re¬ 
sults for a broad range of values of the lattice spacing a. Matters 
are different in the pionful theory, which to LO consists of two 
four-nucleon contact interaction and the long-ranged static one- 
pion-exchange potential (OPEP). As will be shown, combining 
the smearing of the contact interactions with a position-space 
regularization of the OPEP will again lead to results largely in¬ 
dependent of a for the physically sensible range of lattice spac¬ 
ing. Hence, one could use this modified leading-order approach 
to improve the current auxiliary field Monte Carlo simulations 
in Nuclear Lattice EFT. In principle, now it is possible to con¬ 
sider the continuum limit a —> 0, however, we refrain from 
doing that here, as it is sufficient to demonstrate lattice spacing 
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independence for a physically relevant range of a. 


from 1 to 3. The corresponding lattice density operators read 


2. The lattice Hamiltonian 

To set the stage and to introduce our notations, we first dis¬ 
cuss the free Hamiltonian. Its discretized form reads 


P a ’“ (n) = ajj (n) a u] (n) , (4) 

pf’“ (n) = atj (n) T Ijr a u (n) , (5) 

Psf (") ~ a lj (”) vSji'TijfCiu (n), (6) 


n,i,j i 


Hf ™ = i 2 21 2oj ° a l in) a u (n)+ Z (_1) " 


k=\ 


(1) 


x <jj k [alj («) a l } (/I + kl) + a]j (n) a tj ( n - £/)]} . 


Here, a. ..at are the fermionic annihilation and creation op- 
erators with spin and isospin indices i,j, respectively, m N - 
(m p + m n )/2 denotes the nucleon mass and / is a unit vector in 
spatial direction. The summation is over all lattice points n on 
the L 3 lattice. We use a stretched (9(«"'j-iniproved action and its 
coefficients io k are summarized in Tab.Q] see e.g. Ref. SEi. 
m indicates the number of hopping points beyond next-neighbor 
interaction used in the Laplacian discretization in each spatial 
direction and we use m = 4 throughout this paper. The stretch¬ 
ing factor N is introduced to minimize the errors arising from 
the discretized dispersion relations on the lattice especially for 
large momenta where the discretization does not approximate 
the continuum relation E — p 2 /(2;;i,v) anymore. While there 
is some arbitrariness on the exact choice of N depending on 
the values of the respective momentum, N = 3.5 is a sensible 
choice. 
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Table 1: Coefficients for the lattice discretization of the Laplacian, the disper¬ 
sion relation and momentum components depending on the stretching factor 
N. 


The interaction potential consists of two/three terms in the 
pionless/pionful theory at lowest order. The contact interaction 
consists of two terms which can be chosen as 


and Gs[S 2 («) represents the pion propagator times the pion- 
nucleon vertex and is defined as 


G St s 2 (») = X 


exp (-//>■ n)v(p Sl )v(p S2 ) 


(7) 


1 + l;ZZ/ _1) * cos(fepz) 


k= 1 / 


with q n = M\ + 6 loo. v(ps,), v(ps 2 ) are the discretized mo¬ 
mentum components of the first and second pion field which 
yields v(p/) = o\ sin(pi) - 02 sin( 2 /?/) = pi( 1 + 0(p 4 )) with the 
coffecients summarized in Tab. Q] We only use a Oia 1 ) dis¬ 
cretization, because we do not want to expand the respective 
interaction too much. A further improved momentum approx¬ 
imation is linked to a further expanded derivative in position 
space including more interactions at distinct lattice points and 
the locality of the pion-nucleon interaction is lost. These mo¬ 
menta arise from the pion field derivative in the pion nucleon 
Lagrangian = SaKZF^N^t ■ (tr • V)tt/V. To arrive at 
Eq. |7}, we note that the pion propagator is derived from the 
discrete action for instantaneous pions which takes the form [jfj] 

s nn O/) =(-y + 3w o) V JTi («) 7 n (n) 

(8) 

+ 'Yj (- 0* 0J k n i (») ?t/ (« + kl) . 

nj,k 


This is reparametrized by n' t (n) = \fq^ni(n). Finally, the new 
pion action reads 


Wi) =^<("KW 


1 

+ - 

Qn 


^ (- 1 )* u k n' (n) n\ (n + kl) 

n,l,k 


(9) 



where the terms are summed over all lattice sites n and the 
isospin index I = 1,2,3. These terms appear in both versions 
of the EFT considered here. In the pionful theory, one has in 
addition the one-pion-exchange potential 


7/ope = --^2 X ^j G SiS 2 ( n l - n l) 

^ S 1,52,/ H\ ,^2 


(3) 


with gA the axial-vector coupling constant and F n the pion de¬ 
cay constant. Si, S 2 are the respective spin indices which run 


and the respective pion propagator reads 


D„ ( P ) = 


1 + — ^ ^(-l/cos (kpi) 


qn k =1 / 


Furthermore, we introduce a Gaussian smearing 


y(p) 


f (P) = ~r exp 
Jo 


with a stretched 0(a 4 ) improved discretization 

3 3 


( 10 ) 


( 11 ) 


Kp) = 2^^(-l)*w J tcos(fcp,) =p 2 [l +<9(f 6 )] , (12) 


k =0 /= 1 


2 




















where the error estimation is valid for N - 0 and the coeffi¬ 
cients given in Tab.Q] /o is necessary for normalization reasons 
and is given by /o = (1 /L 3 ) Y, P exp[-£>v(p)/2], This smearing 
modifies the contact interaction in momentum space according 
to 

Hcont = \ / (P ) [CP "'' 0 (p)p Qt ’ fl ip) 

P 

(13) 

+c I Y J pf’ a (p)pf’ a (p) ■ 

i 

Such a smearing was introduced in Ref. |Q| to reduce the effect 
of high momentum contributions and remove the clustering in¬ 
stability of the contact interaction. As the leading order (LO) 
contribution is iterated to all orders, such a smearing sums up 
some of the higher order corrections. All other higher order 
corrections are then treated in perturbation theory (as long as 
the lattice spacing is sufficiently large). Here, we concentrate 
on the lowest order and leave the discussion of the treatment of 
the higher order effects to a later publication. In any case, the 
smearing of the LO contact terms is a useful tool to improve 
the description of the S-waves in a very efficient way without 
including all higher-order terms in a chiral counting scheme. In 
general, it is not necessary that the smeared contact interactions 
for the 'So- and \S i-channel have the same smearing function, 
f(p). However, this is important for auxiliary-field Monte Carlo 
lattice simulations. Without the same smearing function for the 
two channels, the Monte Carlo simulation would have a far big¬ 
ger problem with sign oscillations. For this reason we consider 
only one smearing function for both channels. The fact that we 
can approximately describe the effective ranges for both chan¬ 
nels using only one smearing function can be viewed as a fea¬ 
ture of the approximate SU(4) Wigner symmetry of the two- 
nucleon interaction. 

In continuum chiral EFT it is necessary to regularize the one- 
pion-exchange potential due to its singularity at very small dis¬ 
tances. We will show that such a singularity also appears for 
small lattice spacings and we will regularize it in position space 
as suggested in Ref. is}] 


/>) = 


1 - exp 



(14) 


where the denominator is motivated by !F 9 {exp[-r 2 /(2fr)]} oc 
e\p(-bq 2 /2) and n is a free parameter. 


In the following we study the two-nucleon system for dif¬ 
ferent lattice spacings. Therefore we use the finite-volume for¬ 
mulas for the binding energy and the effective range expansion 
given by Liischer i ll ll [12] 


exp 

E b (L) - E“ = -24|A| 2 - 

p cot 8 0 (p) = ^S(q), n = {^j 


(-J=2jiEfL) 
2 uL 


(15) 

(16) 


with A the normalization constant of the wave function for large 
distances and p = m p m n /(m p + m n ) the reduced mass. The 


energy eigenvalues of the system are linked to the momentum 
squared by E - p 2 /(2p) and S (rf) is the zeta function. Its ex¬ 
pansion for small q is given by 


5 (77) = — + So + s i/? 1 + s 2 r + spf + ... ( 17 ) 

n 

and the coefficients S', can be found in @]. We do calcuia- 
tions at volumes large enough that Liischer’s finite volume for¬ 
mulas are sufficient for infinite volume extraction with negli¬ 
gible residual finite volume dependence so we can entirely fo¬ 
cus on the remaining dependence on the lattice spacing. For 
the scattering state calculations, the finite volume corrections 
are of size exp(-L/R), where L is the box length and R is the 
range of the interactions DJ. For the bound state calcula¬ 
tions, the finite-volume formula in Eq. (IT5l > captures the leading 
exp(- yj-2 pE p L), but there are also smaller corrections of size 
exp(- V2 ^J-2pEgL) 1121. For S -waves, the effective range ex¬ 
pansion is p cot do ip) ~ -1/flj + (r e /2)p 2 with a s the scattering 
length and r e the effective range. Hence, we match our finite 
volume results to the infinite volume effective range parame¬ 
ters and the deuteron binding energy, a \ Sn = -23.76(1) fm, 
n So = 2.75(5) fm, and E b = -2.224575(9) MeV. 

In the fits, different lattice volumes are used to produce 
enough data points to make the respective fit using Eqs. (1151161) 
and predict n Si as well as cu Sl . Finally, one can compare 
the prediction to the experimental values given by a3 S] = 
5.424(4) fm and n Sl = 1.759(5) fm. In the following we 
will repeat this procedure for lattice spacings between a~ l = 
100 MeV and or 1 = 400 MeV, respectively, that is approxi¬ 
mately between 2 fm and 0.5 fm in the pionless theory and for 
a between 0.3 and 2.0 fm the EFT with pions. 


3. The pionless theory 

Initially, we consider the pionless case which works well for 
very low energies and is described by the effective Hamiltonian 

H = H frce + W cont . (18) 

ft was shown in Ref. Cl that for any lattice spacing the 
non-smeared contact interaction cannot reproduce the effective 
range correctly. Hence, we introduce a smearing according to 
Eq. (fill) . The calculation is performed for N — 0 and N — 3.5 
to estimate the influence of the stretching factor in the improved 
action. As we only have to consider the S-wave projection, we 
use the appropriate linear combinations 

C=i(3C, So + C3 Sl ) C,= '- (c, So + C3,,) . (19) 

The results for two stretching factors N — 0 and N - 3.5 are 
shown in Fig. [I] For N - 3.5, the explicit values of the fitted pa¬ 
rameters and the predictions are listed in Tab.[[| The lattice size 
of L — 34,36,38 is motivated by the corresponding physical lat¬ 
tice length of Lphy S = 16.66,17.64,18.62 fm for a minimal lat¬ 
tice spacing of a - 0.49 fm. This should be still large enough 
so that higher order terms to Eqs. (1151161) and in the effective 
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a~ l [MeV] 

a [fm] 

C 3Sl [MeV] 

C. So [MeV] 

b [MeV' 2 ] 

a 3 Sl [fm] 

n Sl [fm] 

100 

1.97 

-46.66 

-28.18 

3.83 ■ 10- 5 

5.611(1) 

2.029(1) 

130 

1.52 

-50.71 

-30.16 

3.93 • 10~ 5 

5.630(1) 

2.047(1) 

160 

1.23 

-50.26 

-29.84 

4.04 ■ 10- 5 

5.624(1) 

2.034(1) 

200 

0.99 

-49.89 

-29.70 

4.13- 10- 5 

5.614(1) 

2.017(1) 

240 

0.82 

-50.17 

-29.92 

4.19- 10- 5 

5.607(1) 

2.008(1) 

300 

0.66 

-50.63 

-30.24 

4.23 ■ 10- 5 

5.602(1) 

2.005(1) 

350 

0.56 

-50.85 

-30.38 

4.25 ■ 10- 5 

5.601(1) 

2.008(1) 

400 

0.49 

-50.91 

-30.43 

4.27 ■ IQ- 5 

5.605(1) 

2.015(1) 


Table 2: Fit results in the Gaussian-smeared 
pionless EFT case for L = 34,36,38 and 
N = 3.5. The values for a 3$ and n Sl are 
predictions (modulo the consistency condition 
Eq. (13). 


range expansion are negligible. The LECs in the respective ta¬ 
ble and plot are obtained after rescaling C3 Sl n So —> Ci^/i ^/a ’ 
so as to account for the different volumes. 



tf[fm] 

Figure 1: Low-energy coupling constants and the smearing parameter b in the 
pionless EFT for a lattice size of L = 34,36, 38 and two different stretching 
factors. 


First, we see the parameters belonging to N - 0 and N - 3.5 
have similar values and they approach a constant value when 
the lattice spacing is minimized. But there is a preference for 
N = 3.5 case because for this case, the variation of the LECs is 
significantly decreased for various lattice spacings than in the 
N - 0 case reflecting the improved dispersion relation approx¬ 
imation which is necessary as we receive data points in differ¬ 
ent regimes for different lattice spacings. Therefore, we use 
N - 3.5 for the following calculations. 


Comparing the theoretical predictions with the experimental 
values gives the impression of a large deviation particularly for 
r3 S[ . As we only work at lowest order, we cannot expect to get 
the physical value already. But we can perform a consistency 
check for our predictions as there exists a relation between the 
binding energy, the scattering length and the effective range, see 
e.g. Ref. ill 


E b * - 


1 



( 20 ) 


This equation is valid for positive scattering length a s and the 
expansion in r e la s , with r e the effective range, is useful pro¬ 
vided \r e \«. a s . As we do not include any partial wave mix¬ 
ing, we can use the expansion up to 0{r~/a 2 s ) and compare the 
physical binding energy in this scheme with our lattice pre¬ 
dictions. Then the binding energy is Eb = -2.052(1) MeV 
for the physical values of the scattering length and the effec¬ 
tive range and our lattice predictions give a binding energy of 
E b - -2.009(1)... - 1.999(1) MeV showing that the relation 
in Eq. (l20l > is approximately fulfilled. 

The constant value of b can be explained by the smearing 
function itself and the Fourier transformation of the contact in¬ 
teraction given in Eq. (fT3l> in dependence of the cutoff/ lattice 
spacing. While we work on the lattice, the Fourier transfor¬ 
mation is limited by the maximum momentum n/a. As we 
consider the continuum limit, new contributions to the summa¬ 
tion are added in Eq. (fl3l > as n/a —> oo and the value of the 
normalization constant changes. The normalization constant /o 
approaches a finite limit and the new contributions are expo¬ 
nentially suppressed. Hence, the lattice spacing dependence of 
the smearing constant should approach a constant value as long 
as other discretization errors are negligible. This is also a good 
check for other regularizations procedures. We will use this 
condition to find the best regularization scheme in the upcom¬ 
ing section. It is quite remarkable that the pionless theory can 
be regularized by such a constant Gaussian smearing leading to 
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a-independent results for 0.5 < a < 2.0 fm. Given that the typi¬ 
cal extension of a nucleon is given by a scale of r ^ 0.85 fm, this 
means that the EFT can be used for all momenta that do not lead 
to a resolution of the internal nucleon structure, at least in the 
pionless theory. A direct comparison with the results of a con¬ 
tinuum calculation is difficult because the occuring divergences 
are usually treated in a different way. As done in Ref. 10, 
one can calculate the scattering matrix using the bubble chain 
summation with a regularization f(p) - e\p(-bp 2 /2) similar 
to the smeared contact interaction on the lattice instead of di¬ 
mensional regularization or a finite cut-off. Then the T-matrix 
is expanded up to and including 0(p 2 ) and matched to the ef¬ 
fective range expansion parameters. The values one obtains are 
of the same order as the lattice values but do not match them 
exactly. Note further that the extension to the three-particle sys¬ 
tem is not straightforward because of the Efimov effect that re¬ 
quires the inclusion of a three-body contact interaction already 
at leading order 0. 

4. The pionful theory 

Including pions is necessary for an effective field theory at 
higher energies. Therefore pions are included according to 
Eq. CB and the full Hamiltonian is 

H — //f ree + // con | + //()l’l . (21) 

A problem is caused by the singularities which arise in the 
short-range region of the OPE contribution. On the one hand, 
this singularity exists for any lattice spacing as the relative dis¬ 
tance r = 0 is possible and gives a very large but still finite 
contribution. On the other hand, minimizing the lattice spacing 
leads to additional lattice points with non-zero but very small 
distances to the origin and which give an additional large short- 
range contribution. 

Firstly, we switch off the smearing of the contact interaction 
and include the pion-nucleon interaction according to Eq. @. 
In this case, the predictions for the effective range parameter are 
not close to the physical value and it is still necessary to include 
the smearing of the contact interaction |@]. 

In the case of smeared contact interaction and non-smeared 
pion-nucleon interaction the LECs do not give a reasonably 
close value for the effective range in the 3 5 i-channel. The pre¬ 
diction is reasonable and constant for a lattice spacing larger 
than 1 fm, but it does decrease towards zero for smaller lattice 
spacings. As we do not include a regularization, the divergent 
/rN-contribution for small lattice spacings is more and more re¬ 
solved for smaller lattice spacings resulting in a very sharp po¬ 
tential. The contact interaction in combination with the smear¬ 
ing factor b can compensate this effect but as the 1 S q- and 3 S\- 
channel do not have exactly the same dependence, while it is 
still possible to fit to one of the effective ranges, the agreement 
for the other one gets worse. 

Now we turn on the smearing of the pion-nucleon interaction. 
While a smearing in momentum space according to Eq. (fill) 
could be possible, we follow the arguments given in Ref. |8|] 


and introduce the regularization in position space as proposed 
in Eq. (IT4l >. Then, the new OPE potential reads 

2 

Hope = Yj Yj " 1t2> 

n S 1 , 52 ,/ « 1,»2 

xG SiS! («i -/i 2 )p“;“(«i)pf;“(n2) ■ 

( 22 ) 





Figure 2: Lattice space dependence of various parameters. Upper panel: The 
finite-volume deuteron bound state energy. Central panel: The effective range 
expansion for the ^o-channel. Lower panel: The effective range expansion for 
the 3 Si -channel. 

Throughout this paper, we use n — 4 in Eq. i l I -H i. but we also 
checked for different n and the results are similar. Higher val- 
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a~ l [MeV] 

a [fm] 

C 3Sl [MeV] 

C, So [MeV] 

b [MeV" 2 ] 

a 3 Sl [fm] 

n Sl [fm] 

100 

1.97 

-54.07 

-36.11 

2.59 ■ 10- 5 

5.470(1) 

1.818(1) 

130 

1.52 

-67.11 

-45.47 

2.40 ■ 10- 5 

5.513(1) 

1.878(1) 

160 

1.23 

-69.31 

-46.52 

2.41 • 10- 5 

5.527(1) 

1.899(1) 

200 

0.99 

-63.83 

-41.39 

2.65 ■ 10- 5 

5.523(1) 

1.893(1) 

240 

0.82 

-61.55 

-39.23 

2.83 ■ 10- 5 

5.511(1) 

1.876(1) 

300 

0.66 

-62.20 

-39.44 

2.92 ■ 10- 5 

5.498(1) 

1.858(1) 

350 

0.56 

-63.30 

-40.19 

2.92 ■ 10- 5 

5.491(1) 

1.842(1) 

400 

0.49 

-64.31 

-40.97 

2.91 ■ 10~ 5 

5.491(1) 

1.842(1) 

500 

0.39 

-65.47 

-41.90 

2.89 ■ 10~ 5 

5.500(1) 

1.836(1) 

600 

0.33 

-66.14 

-42.71 

2.82 • 10- 5 

5.553(1) 

1.835(1) 

700 

0.28 

-65.61 

-42.64 

2.86 ■ 10- 5 

5.649(1) 

1.845(1) 


Table 3: Fit results in EFT with Gaussian- 
smeared contact interaction and position space 
regularization for the pion-nucleon interaction 
with N = 3.5. The Lattice size is L = 
34,36, 38 for a = 0.49... 1.97 fm and L = 
38,40,42 for a = 0.39... 0.28 fm. The values 
for a3 S] and rs S] are predictions (modulo the 
consistency condition Eq. ED). 


ues of n are only necessary for strongly divergent contributions 
like two-pion-exchange potentials which we neglect in this ex¬ 
ploratory study. 

In Fig. [2] the finite volume binding energy of the deuteron and 
the respective effective range expansion for the S-waves for a 
subset of the lattice spacings used are shown. Finite volume ef¬ 
fects modify the binding energy so that a correction according 
to Eq. < IT5l) is necessary but higher order terms are still negligi¬ 
ble. Also the data points we obtain for the effective range ex¬ 
pansion still have small enough momenta squared so that the ex¬ 
pansion up to 0(p 2 ) is feasible. The resulting LECs, the smear¬ 
ing constant as well as the predictions for L = 34,36,38 are dis¬ 
played in Tab. [3] While the value of the LECs are different to the 
pionless case, their general shape depending on the lattice spac¬ 
ing does not change and the LECs remain negative as well. Also 
the smearing factor b remains in a certain range between 2.40 
to 2.92 x 10 5 MeV 2 , indicating that this regularization scheme 
works quite well. We do not show the results for the smaller set 
of lattice sizes, but we have shown that finite volume effects 
between these two sets become negligible for lattice spacings 
larger than 0.7 fm. Assuming a large enough lattice there is a 
rise in the smearing constant between 1.8 fm and 0.7 fm, and it 
is constant again for lattice spaces smaller than 0.7 fm. While 
in the range from 0.7 fm and 1.8 fm the pion-nucleon contri¬ 
bution changes, the lattice sites closest to the origin contribute 
more and more to the potential due to the divergent structure of 
the potential. This effect is compensated by the regularization 
at a certain range as we further decrease the lattice spacing, 
and then the regularized potential does not change its shape 
anymore. Note that in the intermediate range there are lattice 
artifacts in the OPEP which cause the oscillatory behavior of 
the LECs as they appear at different physical lengths and multi¬ 
plied by a different regularization factor. As a result, we cannot 
see a plateau for the LECs Cs$ and C'i ,y (1 at a lattice spacing 
of a — 0.5 fm, but we have to further decrease the lattice spac¬ 
ing. Simultaneously, we increase the number of lattice points to 
keep the physical lattice size. The respective values are shown 
in Tab.[3]as well and there is a plateau starting from a = 0.39 fm. 
The predictions for the effective range expansion parameters 
are also quite good as a s Sl as 5.470(1).. .5.649(1) fm and 
ris, ~ 1.818(1)... 1.899(1) fm depending on the lattice spac¬ 
ing. Repeating the calculation for the binding energy according 


to Eq. (l20l > up-to-and-including terms of 0(j"/a]) yields for the 
binding energy Eg as 2.026(1).. .2.033(1) MeV for the vari¬ 
ous lattice spacings in the range between 0.5 fm and 2.0 fm. 
This is again close to the physical binding energy and the re¬ 
spective relation is fulfilled. For the remaining ones, it yields 
Eg as 2.034(1)... 1.899(1) MeV indicating that larger volumes 
become necessary particularly for the smallest lattice spacing. 
Finally, the theory is well regularized and the differences be¬ 
tween the lattice predictions and the experimental values of 
a 3 Sl and r> v , respectively, can be compensated by higher or¬ 
der terms. 

5. Conclusion and Outlook 

In this paper, we have studied the lattice spacing a depen¬ 
dence of the neutron-proton system at leading order in the pio¬ 
nless and the pionful theory. We have used the scattering length 
and the effective range in the 1 S o partial wave together with the 
deuteron binding energy to fix the two LECs related to the LO 
four-nucleon contact interaction and the smearing parameter b. 
To focus on the lattice spacing dependence, we have worked 
at sufficiently large lattice volumes so that finite volume ef¬ 
fects do not play any role here. In the pionless case, it is suf¬ 
ficient to smear the contact interactions with a Gaussian-type 
function, cf. Eq. (ITTI) . to achieve a-independence in the range 
0-5 1$ a < 2.0 fm. We have explicitly shown this for the scat¬ 
tering length and the effective range in the ] S [-channel, being 
aware of the strong correlation between Eg and a , ■ In the 
pionful theory, the contribution from the one-pion-exchange is 
best to be regularized in position space, as discussed in detail 
in Ref. @]. Again, we can demonstrate lattice spacing indepen¬ 
dence for the same range of a. Therefore, it should be possible 
to calculate the phase shifts with the transfer matrix method 
and the spherical wall method. Smaller lattice spacings should 
minimize the lattice errors such as the broken rotational invari¬ 
ance and make it possible to increase the accuracy especially 
for higher momenta in the partial wave decomposition. Clearly, 
when decreasing the lattice spacing, one has to be aware that 
the perturbative treatment of the NLO and higher order correc¬ 
tions becomes doubtful and requires a separate investigation but 
it does not invalidate the results found here. Furhermore, one 
should implement this regularization in Monte Carlo simula- 
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tions in order to minimize the lattice spacing and artefacts in 
the simulation. This, however, is a separate issue. 
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